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ABSTRACT 

Evidence suggests that gamma-ray burst (GRB) ejecta are likely magnetized, although the degree of 
magnetization is unknown. When such magnetized ejecta are decelerated by the ambient medium, the 
characteristics of the reverse shock emission are strongly influenced by the degree of magnetization. 
We derive a rigorous analytical solution for the relativistic 90° shocks under the ideal MHD condition. 
The solution is reduced to the Blandford-McKee hydrodynamical solution when the magnetization 
parameter a approaches zero, and to the Kennel-Coroniti solution (which depends on cr only) when 
the shock upstream and downstream are ultra-relativistic with each other. Our generalized solution 
can be used to treat the more general cases, e.g. when the shock upstream and downstream are 
■ mildly relativistic with each other. We find that the suppression factor of the shock in the strong 

magnetic field regime is only mild as long as the shock upstream is relativistic with respect to the 
downstream, and it saturates in the high-cr regime. This indicates that generally strong relativistic 
shocks still exist in the high-cr limit. This can effectively convert kinetic energy into heat. The overall 
efficiency of converting ejecta energy into heat, however, decreases with increasing cr, mainly because 
the fraction of the kinetic energy in the total energy decreases. We use the theory to study the reverse 
shock emission properties of arbitrarily magnetized ejecta in the GRB problem assuming a constant 
density of the circumburst medium. We study the shell-medium interaction in detail and categorize 
various critical radii for shell evolution. With typical GRB parameters, a reverse shock exists when 
cr is less than a few tens or a few hundreds. The shell evolution can be still categorized into the thick 
and thin shell regimes, but the separation between the two regimes now depends on the a parameter 
and the thick shell regime greatly shrinks at high-cr. The thin shell regime can be also categorized 
' into two sub-regions depending on whether the shell starts to spread during the first shock crossing. 

Q-f The early optical afterglow lightcurves are calculated for GRBs with a wide range of a value, with 

the main focus on the reverse shock component. We find that as a increases from below the reverse 
shock emission level increases steadily until reaching a peak at cr<, 1, then it decreases steadily when 
cr > 1. At large cr values, the reverse shock peak is broadened in the thin shell regime because of 
the separation of the shock crossing radius and the deceleration radius. This novel feature can be 
regarded as a signature of high cr. The early afterglow data of GRB 990123 and GRB 021211 could be 
understood within the theoretical framework developed in this paper, with the inferred cr value >0.1. 
The case of GRB 021004 and GRB 030418 may be also interpreted with higher cr values, although 
more detailed modeling is needed. Early tight optical upper limits could be interpreted as very high 
cr cases, in which a reverse shock does not exist or very weak. Our model predictions could be further 
tested against future abundant early afterglow data collected by the Swift UV-optical telescope, so 
that the magnetic content of GRB fireballs can be diagnosed. 

Subject headings: gamma-rays: bursts — radiation mechanisms: non-thermal — shock waves — stars: 
magnetic fields 
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1. INTRODUCTION 

Extensive broad-band observational campaigns and 
theoretical modeling of gamma-ray burst (GRB) after- 
glows have greatly advanced our understanding of these 
mysterious cosmic explosions. Yet, the origin of the GRB 
prompt emission itself and the nature of the relativis- 
tic flow (which are directly connected to the function of 
the central engine) is still unknown (e.g. Meszaros 2002; 
Zhang & Meszaros 2004). In particular, it is unclear 
how important the role of magnetic fields is in producing 
GRBs. Recently two independent pieces of evidence sug- 
gest that the GRB central engine is likely strongly mag- 
netized. First, the claimed detection of the very high 
degree of linear polarization of gamma-ray emission in 



GRB 021206 (Coburn & Boggs 2003, see however Rut- 
ledge & Fox 2004), if true, could be readily interpreted 
by assuming that the magnetic field involved in the syn- 
chrotron radiation is globally ordered (e.g. Waxman 
2003; Granot 2003), although some alternative explana- 
tions remain (e.g. Waxman 2003). Second, recently we 
(Zhang, Kobayashi & Meszaros 2003, hereafter ZKM03) 
developed a method to perform combined reverse and 
forward shock emission study for GRB early optical af- 
terglows, and revealed that a stronger magnetic field in 
the reverse shock region than in the forward shock region 
is needed to interpret the early afterglow data of GRB 
990123 and GRB 021211. This claim was confirmed by 
independent detailed case studies for both bursts (Fan 
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et al. 2002; Kumar & Panaitescu 2003). These find- 
ings suggest that magnetic fields may play a significant 
role in the GRB physics, as has been suggested by vari- 
ous authors previously (e.g. Usov 1994; Thompson 1994; 
Meszaros & Rees 1997b; Wheeler et al. 2000; Spruit, 
Daigne & Drenkhahn 2001; Blandford 2002). Within 
the framework of the currently favored collapsar progen- 
itor model for GRBs (MacFadyen & Woosley 1999), the 
ejecta are found to be magnetized when MHD simula- 
tions are performed (Proga et al. 2003). 

The degree of magnetization of the ejecta, however, is 
unknown. This is usually quantified by the parameter a 
(see eq.0 for a precise definition), the ratio of the electro- 
magnetic energy flux to the kinetic energy flux. Current 
GRB models are focused on two extreme regimes. In 
the first regime, it is essentially assumed that the GRB 
fireball is purely hydrodynamical. Magnetic fields are 
introduced only through an equipartition parameter cb 
(which is of the order of 0.001-0.1) for the purpose of 
calculating synchrotron radiation. This is the er — ► 
regime. In this picture, the GRB prompt emission is pro- 
duced from internal shocks (Rees & Meszaros 1994) or 
sometimes from external shocks (Meszaros & Rees 1993; 
Dermer & Mitman 1999). This is currently the standard 
scenario of GRB emission. The second is the a — ► oo 
regime. This is the regime where a Poynting-flux dom- 
inates the flow, and GRB prompt emission is envisaged 
to be due to some less familiar magnetic dissipation pro- 
cesses (e.g. Usov 1994; Spruit et al. 2001; Blandford 
2002; Lyutikov & Blandford 2003). In principle, a GRB 
event could include both a "hot component" as invoked 
in the a — model (e.g. due to neutrino annihilation) 
and a "cold component" as invoked in the a — oo model, 
the interplay between both components may result in a a 
value varying in a wide range (Zhang & Meszaros 2002). 
It is an important but difficult task to pin down the de- 
gree of magnetization of GRB ejecta. 

GRB early afterglow data (especially in the optical 
band) potentially contain essential information to diag- 
nose the magnetic content of the fireball. The reason 
is that an early optical afterglow lightcurve is believed 
to include contributions from both the forward shock 
(which propagates into the ambient medium) and the 
reverse shock (which propagates into the ejecta itself). 
Since the magnetization degree of the ejecta influences 
the emission level of the reverse shock (or maybe even 
the level of the forward shock), by studying the interplay 
between the reverse shock and the forward shock emis- 
sion components, one could potentially infer the degree of 
magnetization of the ejecta. In all the current analyses, 
the reverse shock emission is treated purely hydrodynam- 
ically (e.g. Meszaros & Rees 1997a; Sari & Piran 1999; 
Kobayashi 2000; Kobayashi & Zhang 2003a; ZKM03). 
When confronted with the available early afterglow data 
(four cases so far: GRB 990123, Akerlof et al. 1999; 
GRB 021004, Fox et al. 2003a; GRB 021211, Fox et al. 
2003b, Li et al. 2003a; and GRB 030418, Rykoff et al. 
2004), the model works reasonably well for two of them 
(GRB 990123 and GRB 021211), although a good fit re- 
quires that the magnetic field in the reverse shock region 
is much stronger than that in the forward shock region 
(ZKM03). For the other two, the lightcurves are not 
easy to explain with the simplest reverse shock model. 
On the other hand, GRB ejecta could in principle have 



an arbitrary er value. When a is large, the conventional 
hydrodynamical treatment is no longer a good approxi- 
mation, and a full treatment involving MHD shock jump 
condition is desirable. 

It is generally believed that a GRB involves a rapidly- 
rotating central engine. If the magnetic dissipation pro- 
cesses are not significant, field lines are essentially frozen 
in the expanding shells. The radial component of the 
magnetic field decays with radius as oc i? -2 , while the 
toroidal magnetic field decays as ~ R~ l . At the external 
shock radius, magnetic field lines are essentially frozen 
in the plane perpendicular to the moving direction. The 
MHD shock Rankine-Hugoniot relations are greatly sim- 
plified in such a 90° shock. Such relations have been stud- 
ied extensively before both analytically and numerically. 
Kennel & Coroniti (1984) derived some simplified analyt- 
ical expressions applicable for strong 90° degree shocks 
whose upstream and downstream are ultra-relativistic 
with each other. The model was used to treat the pulsar- 
wind nebula problem. In this regime, the strength of the 
shock is essentially characterized by only one parameter, 
i.e. the a parameter. The conclusion was confirmed later 
by numerical simulations (e.g. Gallant et al. 1992) 

Within the context of GRBs, since a GRB invokes a 
transient release of energy, the ejecta shell has a finite 
width (in contrast to the long-standing pulsar wind). 
Under some conditions, the reverse shock upstream and 
downstream could never become relativistic with each 
other when the reverse shock crosses the ejecta shell. In 
the a = limit, whether the reverse shock becomes rel- 
ativistic depends on the comparison between the time 
scale (T) of the central engine activity (essentially the 
duration of the burst) and the time scale (i 7 ) when the 
mass of the ambient medium collected by the fireball 
reaches l/jo times the mass of the ejecta (e.g. Sari & 
Piran 1995; Kobayashi, Piran & Sari 1999). Both times 
are measured by the observer. The case of T > t 7 is 
called the thick shell regime, and the reverse shock is 
relativistic. In many cases, however, one has T < t 7 . 
i.e. the thin shell regime 1 . The reverse shock is initially 
non-relativistic and only becomes mildly relativistic as 
the shock crosses the shell. For magnetized ejecta (e.g. 
a shell with a finite width but an arbitrary a value), the 
separation between the thick and thin shell regimes be- 
comes more complicated, but the non-relativistic reverse 
shock case (for the thin shell regime) is even more com- 
mon (see ^3.3(1 . The Kennel-Coroniti approximation can 
not be directly used. The theory developed in this paper 
becomes essential to discuss the reverse shock physics in 
this parameter regime. 

In this paper, we present a detailed treatment of re- 
verse shock emission for an arbitrarily magnetized ejecta 
under the ideal MHD condition. The reverse shock emis- 
sion in the mildly magnetized regime is also discussed 
by Fan, Wei & Wang (2004a) recently. Here we will de- 
velop a theoretical framework to include discussions for 
the reverse shock emission in a wider a range, as well as 
for various ejecta-medium interaction parameter regimes. 
We first (§2) present a rigorous analytical solution for 

1 If a GRB contains several well-separated emission episodes, 
the whole burst may be even separated into several discrete shells. 
In such cases, even a long-duration burst may be treated as the 
superposition of several thin shells rather than one single thick 
shell. See Zhang & Meszaros 2004 for more discussions. 
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the MHD 90° shock jump conditions, which is applica- 
ble for an arbitrary a value and for an arbitrary Lorentz 
factor 721 between the upstream and the downstream. 
The detailed derivation and the relevant equations are 
presented in the Appendix A. We then (§3) discuss the 
ejecta - medium interaction within the context of GRB 
fireball deceleration and re- investigate the critical fireball 
radii and re-categorize the thick vs. thin shell regimes. 
This leads to a more complicated picture than in the 
pure hydrodynamical case (Sari & Piran 1995). In §4, 
we calculate the synchrotron emission from the shocks 
under the conventional assumptions about the particle 
acceleration in collisionless shocks, and present the pre- 
dicted GRB early optical lightcurves for a wide range of 
a value. We discuss how the early afterglow data may 
be used to diagnose the magnetic content of GRB ejecta. 
Our results are summarized in §5 with some discussions. 

2. ANALYTICAL SOLUTION OF THE RELATIVISTIC 90° 
SHOCKS 

We now consider a relativistic shock that propagates 
into a magnetized medium. In the following analysis, the 
unshocked region (upstream) is denoted as region 1, the 
shocked region (downstream) is denoted as region 2, and 
the shock itself is denoted as "s" 2 . Hereafter Qij denotes 
the value of the quantity Q in the region V in the rest 
frame of "j" , and Qi denotes the value of the quantity 
Q in the region "i" in its own rest frame. For example, 
712 is the relative Lorentz factor between the regions 1 
and 2, f3\ s is the relative velocity (in unit of the speed of 
light c) between the region 1 and the shock, £?2s is the 
magnetic field strength of the region 2 in the rest frame 
of the shock, while B\ is the comoving magnetic field 
strength in the region 1, etc. The relativistic 90° shock 
Rankine-Hugoniot relations could be written as (Kennel 
& Coroniti 1984) 
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where /3 denotes the dimensionless velocity, 7 = (1 — 
^32^-1/2 dgnotgg the Lorentz factor, and u = f3~f denotes 
the radial four velocity, so that 7 2 = 1 + u 2 . Hereafter, 
n, e, p — (F — l)e denote the number density, internal 
energy and thermal pressure, respectively, and F is the 



adiabatic index. The enthalpy is nm p c 
specific enthalpy can be written as 

f 
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where m p is the proton mass and c is the speed of light. 
It is convenient to define a parameter 
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B 
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2 Notice that such a notation system is only valid for §?? and the 
Appendix A. When discussing the GRB problem, i.e. the ejecta- 
medium interaction (§??), we introduce different meanings for the 
subscript numbers. 



to denote the degree of magnetization in each region. The 
magnetization parameter in the upstream region (eri) is 
a more fundamental parameter, since it characterizes the 
magnetization of the flow itself. We therefore define 3 

a^a x = Bl 2 . (7) 
47roiMi7f s 

In our problem, we are interested in a "cold" upstream 
flow, i.e., e% = pi = 0, so that \i\ = m p c 2 . This is 
the only assumption made in the derivation. After some 
algebra (see Appendix), we can finally write 
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n 2 m p c z ' [ 2iti s (721, cf)U2s (721, cr) 

Here «2s (721,0) is a function of 721 and a, and can be 
solved once 721 and a are known. After some analytical 
treatments of the relativistic Rankine-Hugoniot relations 
feq. |H4| ). one can come up with an equation to solve 
u 2s - Since it is com plicated, we only present it in the 
Appendix fea. |A16|). Once u 2s is solved, we can also 
solve Mis (using eq. |A8j ). i.e., 

"ls(721,Cr) = U 2s (721,CT)721 + [u2 s (721,er) + l] 1/2 (72i-l) 1/2 

The compressive ratio can be derived directly from 
cq.lQJ, i.e. 

K(72i^) + l] 1/2 ,., 2 ,u/j 

7 ; (721-F) 

u 2s (72i,o-) 

(10) 

The main point here is that both e2/n 2 m p c 2 and 712/ni 
can be determined by two unknown parameters, i.e. 721 
and cr, so that when they are given arbitrarily, the whole 
problem is solved. 

In the downstream region, the total pressure includes 
the contribution from the comoving thermal pressure 
p 2 = (f — l)e2 and the comoving magnetic pressure 
Pb.2 = -Bf/87r. The ratio between the magnetic pres- 
sure to the thermal pressure is also a function of a and 
721: 
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where eqs.@ and Q) have been used. 

The correctness of the solution (eq. 
in two asymptotic regimes. 



(11) 



16] ') are verified 



2.1. The cr = regime 

When a = 0, the equation t o so lve u 2s (721,0) 
fea. |A16j ^ is greatly simplified fea. |A25| ). All quanti- 
ties can be expressed as a function of 721 . The solutions 
are 



l 2s ' 



(72i-l)(f-l) 2 
f (2 - f X721 - 1) + 2 



(12) 



3 Notice that the definition is slightly different from that in Ken- 
nel &; Coroniti (1984). We find that our definition allows the pa- 
rameters to be coasted into an analytical form as a function of a 
and 721. 
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The equations I|14I16[I are just equations (3-5) of Bland- 
ford & McKee (1976), and they are the starting point for 
the hydrodynamical analysis of the reverse shock emis- 
sion (e.g. Sari & Piran 1995). Under the limit of 721 ^S> 1 
and f — 4/3 (i.e. the downstream fluid is relativistic), 
the equations l|15f) and l|16(l are reduced to the familiar 
forms of n 2 /ni = 47 2 i + 3, 7i s ~ ^721 and j 2s — 3^/4 
(or u 2s ^ \/2/4). 

2.2. The 721 3> 1 regime 

I n the 721 — > 00 limit, the equat ion for w| s (72i, tr) 
fca. |A16| ') is also simplified ('ea. |A26j % ). The solution of 
u 2s is a function of cr only, which reads 
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Notice that there are two solutions with the term ±V X 
in the numerator of ea. (|17f) . but the minus term leads to 
negative pressure and is un-physical. For a relativistic 
downstream region, i.e. T — 4/3, the solution is reduced 
to 



8cr 2 + IOct + 1 + ^M<7 2 {a + l) 2 + 20cr(o- + 1) + 1 
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Fig. 1. — The variations of six parameters, i.e., U2 S , 7is/72l, 
e2/n2m p c 2 , 712/ni, p b 2/P2, and F, as a function of 721- The thick 
solid line indicates the case for a = 0, which is the Blandford- 
McKee (1976) solution. For Pb,2/P2 (panel e), the a = line is at 
negative infinity. The dashed lines, starting from the one closest 
to the thick line, are for a = 0.01, 0.1, 1, 10, 100, 1000, respectively. 
The parameter e2/n2m p c 2 (random Lorentz factor in the shocked, 
downstream region) is normalized to (721 — 1), and the parameter 
n%/ni (compressive ratio) is normalized to (4721 +3), both are the 
values expected in the a = case. 




This is the eq.(4.11) of Kennel & Coroniti (1984). With 
U2s, one can derive u\ s using equation (5). which depends 
on 721 as well. The quantities e2/n2TO p c 2 and 712/ni can 
be also derived accordingly. In the a = limit, equation 
<j*H3|) is reduced to U2 S — which is consistent with 

the asymptotic results in M2.ll 

2.3. The general cases 

For more general cases with arbitrary values of 721 and 
er, u 2 s (72 i, cr ) has to be solved rigorously. The equation 
■cg.IA Hil i is solved numerically, and the solutions indeed 
show deviations from both asymptotic regimes for arbi- 
trary 721 and a values. 

Figure 1 shows the variations of six parameters, i.e., 
U2s, 7is/72i, e 2 /n 2 m p c 2 , n 2 /ni, Pb, 2 /p 2 , and F (see u 2s and 7i s /72i increase with cr, while both e 2 /n 2 m p 
definition in eq.|33j) as a function of 721. The thick and n 2 /ni decrease with a. For 713/721, as long as 721 is 



Fig. 2. — The variations of six parameters, i.e., U2 S , 7is/72l, 
e2/n2iripC 2 , ri2/ni, Pb,2/P2, and -F, as a function of a. The thick 
solid line is the Kennel-Coroniti (1984) solution, denoting a 721 3> 
1 regime. The dashed lines, starting from the one closest to the 
thick line, are for 721 = 1000, 100, 10, 5, 3, 1.5, respectively. Again 
the parameters e2/n2m p c 2 and n2/ni are normalized to (721 — 1) 
and (4721 +3), respectively. 



solid line indicates the case for a = 0, which is the 
strict Blandford-McKee (1976) regime. The dashed lines, 
starting from the one closest to the thick line, are for 
cr = 0.01,0.1,1,10,100,1000, respectively. In order to 
find out the correction factors to the pure hydrodynam- 
ical case, we normalize e 2 /n 2 m p c 2 and n 2 jn\ with re- 
spect to the cr = case. T = 4/3 has been adopted. We 
find that all the parameters achieve asymptotic values 
when 721 3> 1, and that the asymptotic value depends 
on the value of a. In Figure 2 we plot the variations of 
the same six parameters as a function of cr. The thick 
solid line is the 721 3> 1 Kennel-Coroniti (1984) limit, 
and the other dashed lines, starting from the closest to 
the thick line, correspond to 721 = 1000, 100, 10, 5, 3, 1.5, 
respectively. For 721 > 100, the Kennel-Coroniti approx- 
imation is good enough. 

An obvious conclusion from Figs. 1 and 2 is that all 
the (normalized) parameters are insensitive to 721 (espe- 
cially when 721 > a few), but are sensitive to cr. Both 
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mildly large (e.g. > 3), the ratio is essentially the func- 
tion of a only. It starts from the conventional value v2 in 
the cr ~ regime and increases quickly as cr approaches 
unity, which means that the shock leads the fluid sub- 
stantially (in the upstream rest frame) in the high a 
regime. Both e2/ri2m p c 2 and nijn\ are suppressed when 
a increases, but the suppression factor (with respect to 
the cr = limit) is not large, especially when 721 is not 
too small. For example, for 721 > 3, the suppression 
factor for eijn^ is larger than 0.6, while that for nijn\ 
is larger than 0.4. Furthermore, the suppression factor 
reaches an asymptotic value when a approaches several. 
This result is very interesting, since conventionally it is 
believed that the shock is completely suppressed when 
cr reaches larger values. Our results suggest that rela- 
tivistic strong shocks still exist in the high-a regime. The 
suppression factor, which essentially does not depend on 
cr, is only mild as long as the shock is relativistic. The 
overall efficiency of converting the total energy (the ki- 
netic energy plus the Poynting flux energy) to heat still 
decreases steadily with increasing cr. The reason is not 
that the shock (which converts the kinetic energy into 
heat) itself is less strong, but is that the fraction of the 
kinetic energy in the total energy, i.e., (l+cr) -1 , becomes 
smaller as a becomes larger. 

3. EJECTA - MEDIUM INTERACTION 

3.1. Basic equations 

Now we consider an arbitrarily magnetized flow with 
magnetization parameter a and Lorentz factor 7 = 74 
being decelerated by an ambient medium with density 
n = m. A pair of shocks form when the shock form- 
ing condition is satisfied, i.e. the relative velocity be- 
tween the two colliders exceeds the sound velocity in the 
medium and the magnetoacoustic wave velocity of the 
ejecta, and that the pressure in the shocked region ex- 
ceeds the pressure in the unshocked region. In the GRB 
case (a relativistic ejecta), a forward shock always form, 
while a reverse shock may not always form if a is too 
large. In the high-cr regime, the magnetoacoustic wave 
velocity is essentially the Alfven velocity for a 90° shock. 
The first condi tion fo r the reverse shock formation is 
741 ^ 1A ~ vl + o~, where ja is the Alfven Lorentz 
factor in the ejecta. For GRBs we have 741 > 100 (to 
ensure that the observed gamma-ray spectrum is non- 
thermal), so that this condition is satisfied as long as 
cr < 10 4 . The second condition is generally more strin- 
gent, which is expressed in eqs. I|31|l and (|43|l below. 
In any case, with reasonable parameters a reverse shock 
could be formed when cr is less than hundreds or tens. 
When the reverse shock forms, we can then investigate a 
picture where two shocks and one contact discontinuity 
separate the ejecta and medium into four regions. Below 
we take the usual convention to define the four regions: 
(1) unshocked medium; (2) shocked medium; (3) shocked 
ejecta; (4) unshocked ejecta. Notice that hereafter the 
numerical subscripts have different meanings from the 
ones used in §?? and Appendix, where the numbers "1" 
and "2" denote upstream and downstream, respectively. 
For the forward shock, both sets of notations coincident, 
while for the reverse shock, the previous "1" and "2" are 
replaced by "4" and "3", respectively. We assume that 
the ambient medium is not magnetized so that <j\ = 0, 



but will assign an arbitrary magnetization parameter 

(20) 



cr = (T4 



B 4s 



for the ejecta. Since we are discussing the problem in the 
rest frame of the medium, we will drop out the subscript 
"1" whenever it means "in the rest frame of the region 
1" . We can then write the following relations based on 
the shock jump conditions. Throughout the following 
discussions, T = 4/3 is adopted. 

— ^ = ( 72 -l) ~ 72 , (21) 
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are the correction factors for /n^m p c 2 and n% / 714 with 
respect to the a = limit, and U3 S is a function of 734 an d 
cr, whose solution is found in the Appendix A fea. |A16j 'l. 
The functions f a and fb are plotted as a function of a in 
Fig.J2t) and Fig.J2i), respectively, which shows f a — > 1 
and ft, — > 1 when a — > 0. Constant speed across the 
contact discontinuity gives 

72 = 73 , (27) 

and constant pressure across the contact discontinuity 
gives 

p^ p^ R? Pn r>t n 

(28) 



3 3 8vr 3 { p 3 ' ' 



or 



where 



62 = e 3 / c , 



fc = / c (cr, 734) = 1 + —, (oc cr) 

P3 



(29) 



(30) 



The pressure ratio Pb,3/P3 is calculated according to 
ea. lfTT|l . whose dependence on a is plotted in Fig. (05 ) 
which shows a oc cr dependence in the cr 3> 1 asymptotic 
regime. Hereafter we denote the asymptotic behavior in 
the cr 3> 1 limit in a pair of parenthesis immediately 
following an equation. 

In order to have equation 129fl satisfied, the condition 
should be that the thermal pressure generated in the for- 
ward shock is stronger than the magnetic pressure in the 
region 4. This gives B|/87r < (4/3)7f rtim p c 2 . Notic- 
ing equation (|20() . the condition for the existence of the 
reverse shock can be written as 



8 9 ni 
3 n 4 



(31) 
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When the reverse shock exists, using equations i|2H24[l 
and (|29|) . one can finally get 



"4 F _ (72- 1) (4 7 2 + 3) 
n\ (734 - 1)(4t34 + 3) 



(32) 



where 

F = fafbfc , (33) 

and f a , f and f c are defined in equations (|25|l . (|26|l and 
H30|l . respectively. The parameter F has been calculated 
for different input parameters, and the results are shown 
in Fig. (If) and Fig.(2f), respectively. We can see that F 
is very insensitive to 734, and is essentially a function of 
cr only. The asymptotic behavior in the a 3> 1 regime is 
F(a) oc a. We then write 

F(^a)~F(a),(<xa). (34) 

The equation (|32JI can be used to define whether the 
reverse shock upstream and downstream are relativistic 
with each other. Similar to the analysis of Sari & Piran 
(1995), we analyze the value of (n^/nijF/^ 2 . The rel- 
ative Lorentz factor of the reverse shock upstream and 
downstream is 



734 



74 
73 



73 
74 



734 > 
„4 /„, 4 



(35) 



1 , we have 



For the relativistic case, i.e. 

{ ni /n x )Fhl ~ 7 2 7(73 2 474 2 ).~ 7 2 4 /7 4 4 « 1- On the 
other hand, for a non-relativistic case, we have 734 ~ 1, 
74 ~ 73 and (734 — 1)(4734 + 3) = e <C 1. This gives 
{ri4,/nx)F/"i1 ~ 1/e > 1. We thus conclude that the 
reverse shock upstream is relativistic with respect to 
the downstream when 7! ^ (n^/n^F, while it is non- 
relativistic when 74 <C (n^/n^F. For a = 0, we have 
F = 1, and the result is fully consistent with Sari & 
Piran (1995). 

3.2. Critical radii 

We consider an isotropic fireball with total energy 
E = Ek + Ep, where Ek is the kinetic energy and Ep is 
the Poynting flux energy. The discussions are also valid 
for a collimated jet by regarding the various energy com- 
ponents as the "isotropic" values. With the definition of 
a (eq.jZ|), we find E P /E K ~ cr (from eq.g]) 4 , so that 
E = E K (l + cr) or E K = E/(l + a) (sec also Zhang & 
Meszaros 2002). We follow the traditional convention to 
define the Sedov length I ~ (E /nirripC 2 ) 1 ^ 3 , where the 
total energy is adopted. The shell baryon number den- 
sity 714 is however defined by Ek- The density ratio is 
n 4 /ni ~ Z 3 /[7f Ai? 2 (l + cr)], where A = max(A 0) -R/74) 
is the thickness of the shell, Ao = cT is the initial width 
of the shell, and R is the fireball radius. This holds for 
both a non-spreading shell (where A = Ao is a constant) 
and for a spreading shell (where A ~ R/-f 2 ). 

Below we revisit the four critical radii related to the 
reverse shock deceleration (Sari & Piran 1995). In our 

4 The full presentation of eq.|l] should be £ti«i s + pi/ni«i s + 
B'( s /8nniui s + Bf a f3f a /8nniui a = fi2U2s + P2/ri2U2 a + 
B% s /8nri2U2s + -B| s /3| s /87rn 2 «2 a - So strictly speaking, E P /E K = 
cr(l + /3?)/2/3^. In the high-tr regime, we have f3± s ~ 1 and the 
factor (1 + Pf s )/2f3f s ~ 1. In the low cr regime, /3 ls < 1, but the 
factor [l + cr(l + /32 a )/2/32 s ] is in any case ~ 1 . We therefore neglect 
the (1 + f3f g )/2/3f s factor in the following discussion. 



following discussion, we assume that a reverse shock ex- 
ists. The asymptotic behaviors at a ^> 1 for various cor- 
rection factors (presented in parentheses) therefore arc 
valid for the cr range in which the reverse shock forming 
condition is satisfied. 

1. The fireball radius for the relative Lorentz factor 
between the reverse shock upstream and downstream 
(i.e. 734) to transform from the Newtonian regime to 
the relativistic regime can be estimated according to 
7! ~ (ni/n^F, which gives 



Rn 



t 3/2 



"1Z 



N 



where 

Cn (734, cr) 



C N (a) = 



F(cr) 



1 



1/2 



(36) 



l,(occr°) (37) 



is the correction factor of Rn with respect to the cr = 
case. Since both F(a) and (1 + cr) have the same asymp- 
totic behavior (oc a) at high cr, the final correction factor 
Cat is always of order unity throughout (see Fig|3J), and 
we will neglect it in the following discussions. 

2. The radius where the reverse shock crosses the shell 
is approximately (Sari & Piran 1995) Ra = [A/(/?4 — 
/?2)] [1 — (74/73X^4/713)]. In the cr = case, the factor 
[1 — (74/73)("4/«3)], which delineates the relative com- 
pression factor due to shock crossing 5 , is a factor rang- 
ing from 1/2 — 6/7, which was neglected for order-of- 
magnitude estimates (Sari & Piran 1995). For arbitrary 
cr values, this parameter is cr-sensitive (through the de- 
pendence of /f,(cr), see eq.j^ni), i.e., becomes <C 1 when 
cr ^> 1, so we can not drop it out. Following the similar 
procedure as in Sari & Piran (1995), and replacing n 4 /ni 
(in the cr = case) by (rii/n^F, one finally has 



ni 



1/2 



Ra - 74A 
where 

Ca (734, <j)^Ca{<?) 



1 _ 1^±] „ A 1 /4/3/4 Ca 
73^3/ 

(38) 



F(cr) 



"I 1/4 



1 



74 n 4 



73 m 



1/2 



1 + crJ 

i-^Y /2 , (occr-VS) (39) 



73 "3 



J • 



is the correction factor of Ra with respect to the a = 
case. This correction factor suggests that the reverse 
shock crosses the shell faster when a becomes larger. 

3. The conventional "deceleration radius" (for the thin 
shell case) is modified in the high-cr regime. According to 
eq.(J21), the ratio of the comoving magnetic fields in region 
4 and 3 is B4/B3 = u^s/u^s. The lab- frame ratio of the 
Poynting flux energy in both regions can be written as 



E 



PA 



74 (-Bf/ n 4) _ M 3s 74 

Ep.3 7 3 (B|/n 3 ) u 4s 7 3 



1 



(40) 



where ~ applies in the cr 3> 1 limit so that Ui s ~ 74 S and 
"3s ~ 73s- Since 74 S /73 S = 74/73 (both have the same 
relation with 734 fea.|35|). ea. l40JI is naturally derived. 

5 If one assumes that after shock crossing a shell with width A 
becomes A', this parameter is simply (A — A')/A. 
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We have calculated this ratio numerically and found that 
in the high-u limit, the difference between this ratio and 
unity is a small quantity comparing a^ 1 . This manifests 
that shocks in the high-cx limit only effectively dissipate 
the kinetic energy (in the baryonic component) in the 
upstream, and the Poynting energy (in the lab frame) 
essentially remains the same. As a result, one should 
define the deceleration radius (for the thin shell case) 
using Ex alone, so that 



I 



1 



C 1 , 



2/3/, . u/3 2/3 T 

74 (l + cr) 1/j 74 



where 



C» = (1+a)- 1 / 3 (oc ff "V3) 



(41) 



(42) 



Notice that the radius i? 7 is still the radius where the 
fireball collects I/74 of fireball rest mass 6 . 

In the above discussion, we already assumed that a 
reverse shock exists. In order to satisfy the condition 
iJSTJ at i? 7 , one can give a more explicit constraint on a, 
which reads 



a<ioof^iy 

V300/ 



T 

10s 



3/2 



E 



10 52 ergs 



-1/2 



(43) 



We can see for typical GRB parameters, a reverse shock 
exists when a is smaller than several tens to several hun- 
dreds. 

It is worth noticing that at the deceleration radius, 
the Poynting energy is not transferred to the ISM yet. 
At the end shock crossing, the magnetic pressure behind 
the contact discontinuity balances the thermal pressure 
in the forward shock crossing. It is not until the reverse 
shock disappears during the deceleration phase when the 
bulk of magnetic energy in the ejecta transfers to the 
forward shock region. During the deceleration, the mag- 
netic fields push the contact discontinuity from behind 
and transfer energy through pdV work. Eventually, the 
total energy will be transferred to the ISM, so that the 
late time afterglow level is still defined by the total en- 
ergy of the fireball. The detailed energy transfer process 
is complicated and will be studied carefully in a future 
work (see Zhang & Kobayashi 2005 for a brief discussion). 
This point is relevant to the calculations of the forward 
shock emission level, and it will be further discussed in 

S3 

In FigOl we numerically plot the functions C/v(c), 
Ca(c) and C 1 (a) for both the mildly relativistic case 
(734 = 1.5) and the extremely relativistic case (734 = 
1000). We can see that Cn is insensitive to both 734 and 
a, and we will treat it as a constant of order unity. The 
correction factor Ca is rather insensitive to 734 and is 
essentially a function of a only. In the a 3> 1 regime, we 

6 This could be derived using energy conservation equation in 
the lab frame before and after the shock crossing(s), i.e. 74(Moc 2 + 
U B ,o) + M lSM c 2 = 7 2 (M c 2 + 7 2 M ISM c 2 + U B ), where M is the 
mass in the original ejecta, Mism is the collected ISM mass as 
the reverse shock crosses the shell, £/b,o is the initial comoving 
magnetic energy, and U& is the comoving magnetic energy after 
shock crossing. This gives (74 — 72)Moc 2 + (i?p,4 — Ep^) = 

(7! - 1)A^ismc 2 5 where E P ^ = 74%, and Ep,3 = 72U B - Ac- 
cording to equation I4UI . the term (Bp 4 — Ep^j drops out from 
the energy conservation equation, so that the equation is effec- 
tively the familiar hydrodynamical one with the total energy being 
E K = E/(l + a). 



_C N (°), T 3,-1.5 



Fig. 3. — The functions Cjv(<t), Ca(<t) and C 7 calculated for 
both 734 = 1.5 and 734 = 1000. 



have Ca oc a 1 I 2 . By definition, C 7 is the function of a 
only. 

4. Finally, the radius where the shell spreads is still 
defined by 

R s ~ 7 4 2 A • (44) 
Taking the convention to define (Sari & Piran 1995) 

£ = (ll^ 2 hT , (45) 
one has the following equation 



where 



Rn _ R y _ £ 1/2 „ _ , 2f? 
4 o 7 Ca 

(VAo) 1/2 (h\ 1/2 



£0 



4/3 
74 



(46) 



(47) 



is the £ value for R < R s (i.e. no spreading occurs), 
where 



Rj 



I 



~t ~ C -v 2 r 8/3 ' 

o 7 7 4 c j^' c 



(48) 



Notice that this notation is exactly the same as in the 
a = case, and the total energy E is used (in defining 
I). This makes £ 7 a constant not depending on a, which 
is convenient for the discussions of the various parameter 
regions in the next section. 

3.3. Parameter regions: thick vs. thin shell regimes 

Equating the four critical radii defines six critical lines 
in the £0 — cr or the T/t^ — a space. This is justified by the 
fact that the spreading regime (which makes £ deviating 

from £0 = (t-y/T) 1 / 2 ) always happens below the critical 
lines, and hence, does not influence the location of the 
critical lines. The six critical lines are 



Rj 


'Ra, 


T 

t 7 " 


Rn ~ 


'Ra, 


T_ ^ 


Rn r 


^ Rji 


T 
1- 


Rn ' 


^ R s , 




R~f ' 


^ R s , 


T_ r 


Ra ' 


^ R s , 


T n 



'(&) 

c A - 4/3 
- c- 2 ~ 

f ~ 1, 
<c 7 ~q 



c 



4/3 



- Q, (oc a 2 / 3 ) 

Q, (oc^ 3 ) 
(« a ) ; 
-V2, (oc a' 1 / 3 ) 
r\ (oc a- 2 / 3 ) 



(49) 
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Fig. 4. — Parameter regimes in the (T/t 7 ) — a space. 



We can see that the first three lines have a same asymp- 
totic behavior at high-cx, and calculations show that they 
essentially coincide with each other. Hereafter we define 

Q{a) ee [C»p 2 ~ C A 4/3 ~ (C 7 /C A ) 4 (oc a 2 / 3 )(50) 

so that Ca ~ Q~ 3 ^ 4 and C 7 
be re-written as 



Q 1/2 , and eq.JUJ can 



T 



(51) 



In principle, changing the order between Rn and R s and 
between 7? 7 and R s does not lead to essential modifica- 
tions of the shock crossing and deceleration physics, so 
that the 4th and 5th lines in ea. l|49|l are not crucial. We 
therefore have two essentially lines that separate three 
physical regimes in the T/t 1 — a space (FigQJ. 

Region (I): the thick shell regime. This is the region 
where T/i 7 > Q is satisfied. In this region, one has 
Rn < R 1 < Ra < R s - The downstream becomes rela- 
tivistic with respect to the upstream before the reverse 
shock crosses the shell (Sari & Piran 1995). The full 
deceleration occurs at the end of shock crossing, i.e. at 
R ~ Ra (Kobayashi et al. 1999). Given a certain ob- 
served central engine activity time T, a total energy E, 
and an ambient density n, one can also define a critical 
Lorentz factor 



7c 



1 



3/8 



(52) 



where some refined coefficients and the cosmological time 
dilation factor are explicitly taken into account in order 
to get the numerical value (here z is the GRB redshift). 
For 7q > 7c, one is in the thick shell regime, while for 
7o < 7ci one is in the thin shell regime. The function 
Q{<j) is plotted as the top curve in FigQ] In the a <C 1 
regime, Q(a) ~ 1. For a 3> 1, we have Q(a) oc cr 2 / 3 , and 
hence, 7 C oc cr 1 / 4 . We can see that given same values for 
the other parameters, the parameter space for the thick 
shell regime is greatly reduced when a is high. More 
bursts are in the thin shell regime. 

Region (II): the non-spreading thin-shell regime. This 
region is defined by Q~ x < (T/t y ) < Q, in which Ra < 
i? 7 < Rn and Ra < R s are satisfied. The common 
features of this regime are that the reverse shock crosses 
the shell (at Ra) before noticeable deceleration occurs 
(at i? 7 , e.g. the Lorentz factor is reduced by roughly a 
factor of 2), that the relative speed between the upstream 



and the downstream never becomes relativistic, and that 
the shell does not spread during the first shock crossing. 
The separation between Ra and R 7 leads to some novel 
features for the reverse shock emission. During the first 
shock crossing the shell is heated so that electrons start 
to emit. However, after the first shock crossing, the shell 
is not decelerated significantly. It is difficult to delineate 
the detailed process at this stage, but a rough picture is 
that higher order shocks may form and bounce back and 
forth between the inner edge of the shell and the contact 
discontinuity. This happens until the shell reaches i? 7 . 
It is likely that the shell remains heated by the multi- 
crossing of shocks and electrons keep radiating at a high 
level for an extended period of time. One then expects 
a broad reverse shock emission peak, which is a novel 
phenomenon in the high-c thin-shell regime. Notice that 
when a is very large, a reverse shock may not form at Ra 
at all. However, since Ra is the smallest in the problem, 
whenever a reverse shock forms, it quickly crosses the 
shell in a radius of Ra , and the above discussion is still 
valid. 

Region (III): the spreading thin-shell region. This is 
defined by (T/t 7 ) < in which R s < R A < i? 7 < 

Rn is satisfied. Since Ra < R~y, again multi-crossing 
of shocks are needed to slow down the ejecta, and the 
downstream never becomes relativistic with respect to 
the upstream. The novel feature in this region compared 
with the region (II) is that the shell starts to spread 
before shock crossing, so that the three radii have the 
relationship 



C^ 3 R 



N 



i? 7 



(Cj/Ca) 2 Ra, 



Q ^^Rn ~ i? 7 ~ Q^^Ra ■ 
We can see that the triple coincidence Rn = i? 7 



(53) 

(54) 
= Ra 



in the thin shell regime (Sari & Piran 1995) is only valid 
when a is small. According to FigQ] this practically 
happens when a < 0.01. 

3.4. Critical times 

We finally derive the shock crossing time t x and the 
deceleration time tdec as measured by the observer. In 
the literature to study the a — regime, t^ ec — t x ~ 
Ra/^2 c ~ max(T, i 7 C 7 ) ~ max(T, f 7 ) has been conven- 
tionally adopted (noticing C 7 = 1 when a — 0). When 
an arbitrary a value is adopted, there are further compli- 
cations to quantify these critical times. First, although 
in the thick shell regime (I) tdec = tx is still valid, in 
the thin shell regimes (II and III) the deceleration radius 
R 7 is larger than the shock crossing radius Ra, so that 
tx < tdec- Second, the time scale RaH\c only describes 
the delay time scales for the emission coming from the 
radius Ra with respect to the emission from the internal 
shock radius, for an infinitely thin shell. A more pre- 
cise description of the reverse shock emission peak time 
should include the thickness of the radiation region. The 
real shock crossing time should correspond to the epoch 
when the emission from the end of the shell reaches the 
observer (see FigGDfor illustration). This gives 

tx~^ + -- (55) 

72 C c 

Here A = max(Ao, R/'y 2 ), so that our definition is valid 
throughout the (T/t 7 ) — a space. In the u « 1 limit, 
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R 



is 



p radius 

A 



Fig. 5. — The space-time diagram of the GRB fireball evolution. 
Solid lines are the world lines of both ends of the shell, and the 
dashed lines are the world lines of light. The lowest dashed line is 
the first light (from the internal shock) that reaches the observer, 
and the other two dashed lines indicate the light emitted from both 
ends of the shell at the shock-crossing radius, Ha- The observed 
shock crossing time is the sum of the delay time R/\/ (2y%c) and 
A/c. To order of magnitude estimate, the compression of the shell 
after shock crossing is not taken into account. 



we always have RaH\c ~ A/c so that to order of mag- 
nitude estimate, one can drop the latter term. In the 
a 3> 1 limit, however, in certain regimes one could have 
R\/l2 c A/c. The correction factor introduced here is 
therefore essential to delineate the reverse shock behav- 
ior in the high-cr regime. We notice that Nakar & Piran 
(2004) recently also noticed this correction within the 
context of a — regime, although their eq.(2) is slightly 
different from our definition. 

In the thick shell regime (I) the reverse shock is 
relativistic at the crossing radius, and one has 7! = 
7 4 [(n 4 /ni)F] 1 / 2 . Using the definition of R A (eq.jSHl), 
one has RaH\c ~ TCl ~ TQ~ 3 / 2 . Since A/c = 
A /c = T, with ca. ljSa t we have i x (I) = t dec (l) ~ 
T(l + <3~ 3 / 2 ). In the region (II), i.e., the non-spreading 
thin shell regime, one has 72 ~ 74 and A = Ao. Using 
eqs.(H)| and gSJ, we get i x (II) - t^T^CA + T ~ 
4 /4 T 1 / 4 g- 3 / 4 + T. In the region (III), i.e. the spreading 
thin shell regime, one has 72 ~ 74 and A/c = RaH\c- 
With eq.JSH, one has f x ~ 2t^(C A /C 7 ) ~ 2i Y Q~ 1 . 

We define the deceleration time tdec as the epoch when 
the fireball is significantly decelerated. For the thick shell 
regime (I), this coincides with the shock crossing time. 
For the thin shell regimes (II, III), it is defined when 
the fireball reaches i? 7 . Following the same argument to 
derive (|55|l (FigJSJ), one can generally define 



tde 




try (.Jry 



(I) 
(II, III) 



(56) 



where we have interchanged 72 and 74 for the thin shell 
regimes. For the regime II, one has t<j e c(H) ~ tyQ -1 / 2 + 
T, while for the regime III, one has idec(IU) ~ t 7 (Q~ 1 / 2 + 

Q- 1 )- 



For convenience, we collect the critical times as the 
following, 

r(i + cr 3 / 2 ), (i) 

t 3 / 4 T 1 / 4 <2~ 3 / 4 + T, (II) , (57) 
2i 7 Q-\ (III) 



and 



r(i + er 3 / 2 ), (i) 

tdec ~ { hQ- 1 ' 2 + T, (II) 

tyiQ-^ + Q- 1 ), (HI) 



(58) 



As demonstrated above, in the thin shell regime, the re- 
verse shock emission should show a broad peak due to 
the separation between t x and tdec- The width of the 
peak of the reverse shock peak is defined as 



{tdec t x ) 



0, (I) 
MQ-V 2 - (T/t 7 ) 1/4 Q _3/4 ], (II) 

t-t(Q~ 1/2 - Q^ 1 ), (HI) 



(59) 



4. SYNCHROTRON EMISSION AND EARLY AFTERGLOW 
LIGHTCURVE 

4.1. Particle acceleration 

The hydrodynamical solution presented above is 
generic and lays a solid foundation for calculating syn- 
chrotron emission in the reverse shock and the early 
afterglow lightcurves. In this section, we will turn to 
the less certain aspects of the problem, i.e. particle 
acceleration and synchrotron emission from the reverse 
shock. In the conventional a = models, particles (ions 
and electrons) are assumed to be accelerated from the 
collision-less relativistic shocks through the first-order 
and probably also second-order stochastic Fermi accel- 
eration mechanisms (Fermi 1949; Blandford & Eichler 
1987). In the standard afterglow models, accelerated 
electrons are assumed to have a power law distribution 
with dN e /dj e oc 7~ p . Numerical simulations of parti- 
cle acceleration in relativistic shocks confirm this simple 
treatment (e.g. Gallant et al. 1992; Achterberg et al. 
2001) in the a = limit. For MHD shocks as discussed 
in this paper, more physical processes enter the prob- 
lem. For example, in the 90° shock problem discussed 
in this paper, the existence of the electrostatic potential 
in the shock front plane and the influence of the Lorentz 
force exerted on the particles by the magnetic and elec- 
tric fields in the upstream tend to trap ions in the shock 
plane. The influence of these effects on shock accelera- 
tion is unclear, although some investigations in this di- 
rection have started (e.g. Double et al. 2004; Spitkovsky 
& Arons 2004) . Detailed treatments of particle acceler- 
ation in the high-cr regime is beyond the scope of the 
current paper. Lacking a detailed model, here we simply 
extend the approach used in the a = regime to arbi- 
trary a values, i.e. we assume a power law distribution of 
electron energy and assign the equipartition parameters 
e e and cb for the electrons and magnetic fields, respec- 
tively. Such an approximation is proven valid when a is 
small, but may progressively become not good enough as 
a increases, especially when a achieves very large values. 
The lightcurves we calculate below nonetheless provide 
a first-order picture on how the reverse shock emission 
level depends on a. 
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4.2. Magnetic fields 

We follow our previous approach (ZKM03) to compare 
the flux level between the reverse shock peak and the 
forward shock peak. This is because, when studying the 
ratio of the peak flux levels, only the ratios of the micro- 
physics parameters (e.g. g(p) = (p - 2)/(p - 1), e e , e B , 
etc.) matter, and one does not need to invoke the ab- 
solute values of those parameters which arc rather un- 
certain. A crucial parameter to study synchrotron spec- 
trum is the comoving magnetic field strength in both 
shocked regions. For the forward shock region, since the 
medium is usually not magnetized, the downstream mag- 
netic field is usually quantified by a fudge parameter ebj, 
which reflects the strength of the magnetic field presum- 
ably generated in-situ due to a certain plasma instability 
(e.g. Medvedev & Loeb 1999). This magnetic field is 
randomly oriented, as has been supported by the ob- 
served weak-polarization level for the optical afterglow 
emission (e.g. Covino et al. 2003 for a review). The 
strength of this magnetic field component is low, with 
€bj ~ (0.01 — 0.001), as inferred from broadband after- 
glow fits (Panaitescu & Kumar 2002; Yost et al. 2003). 
For the reverse shock, in the current model the magnetic 
field is predominantly due to the compression of the up- 
stream magnetic field. Its level depends on the a value of 
the upstream, and the field is globally structured, so that 
the optical flash due to the reverse shock emission (such 
as the ones observed from GRB 990123 and GRB 021211) 
should have been strongly polarized (see also Granot & 
Konigl 2003; Fan et al. 2004a; Sagiv, Waxman & Loeb 
2004). 

The forward shock comoving magnetic energy density 
is defined by 



8^ = 8^ = e2£BJ 



(60) 



where ebj is the conventional magnetic equipartition 
parameter in the afterglow theory which delineates the 
fraction of the total internal energy that is distributed to 
magnetic energy. In the reverse shock region, the comov- 
ing magnetic energy density is dominated by the shock- 
compressed upstream magnetic field, and can be denoted 
as 

l]1 * C/c-1)^, (61) 



8tt 



where f c = 1 + Pb,z/pz (eq. [201) has been used. For easy 
comparison (with respect to the conventional definition 
of e_B,/), we can write Ij61|) as 
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C2tB, 



where 



_ (/c-1) 

3/c 



(62) 



(63) 



is an artificial parameter to simplify the discussions. 
With this definition, we can write 



B f \£B,f 



1/2 



(64) 



This is an important parameter which delineates the ra- 
tio of the magnetic field strength in the reverse shock 
and forward shock regions. This ratio has been found to 




Fig. 6. — The equivalent "magnetic equipartition parameter" 
in the reverse shock region, i~B,r, as a function of cr, calculated for 
both 734 = 1000 and 734 = 1.5. The dashed line is the assumed 
magnetic equipartition parameter in the forward shock region, i.e., 
e B,f ~ 0.001, which is also the "bottom level" for the reverse 
shock magnetic field in the low-cr reg ime. The thin solid lines 
are calculated completely from eg. 1631 . while the thick solid lines 
include the contribution of the random field in the low-cr regime. 
The fact that the parameter eg r approaches 1/3 in the high-cr 
regime i s on ly an artificial effect resulting from the definition of 
e B , r (eq.E3). 



be larger than unity in GRB 990123 and GRB 0201211 
(e.g. ZKM03), and discussion of this parameter is very 
essential to quantify the relative emission properties of 
both shocks. 

In Fig E| we plot 1b, r as a function of a for both 
734 = 1000 and 734 = 1.5. We can see that it in- 
creases with a initially, and saturates at a value 1/3 as 
(j > 1. The asymptotic behavior is already obvious in 
cq. I|63|) . We note that the number 1/3 is a pure artificial 
effect given the definition of eb,t (eq.jfH]). The "real" 
magnetic equipartition factor in the region 3 approaches 
unity when a 3> 1. In the figure we also plotted the 
magnetic equipartition parameter in the forward shock 
region, i.e., esj ~ 0.001. This level could also be re- 
garded as the "bottom level" in the cr -C 1 regime for 
the reverse shock region. The thick lines in Fig|S] arc 
the "total" eb in the reverse shock region (which include 
both the amplified structured field component and the 
random field component), which saturates to es.f in the 
low-cr regime. 

An interesting conclusion drawn from FigO is that 
no matter what a or 734 values are taken, the ratio 
£b/cbj is at most ~ 300 for ebj ~ 0.001. Or effec- 
tively, the reverse-to-forward shock magnetic field ra- 
tio TZb can not be significantly larger than 15. Since 
this value was inferred from the case of GRB 990123 
(ZKM03), we tentatively conclude that the ejecta in 
GRB 990123 has cr > 0.1 where eb,t reaches its peak 
value. The absolute values of eb,t and ebj are consis- 
tent with those obtained from the detailed modeling. For 
example, cbj ~ 7.4 x 10~ 4 was inferred by Panaitescu 
& Kumar (2002), while e B , r /e B j ~ 15 2 ~ 225 was in- 
ferred by ZKM03, so that e B , r ~ 0.17. This is close to 
the maximum es,r we have calculated. 

4.3. Lightcurve peak times and peak fluxes 

Before describing the detailed process of calculating 
afterglow lightcurves, it is informative to define the so- 
called peak times and peak fluxes. An early optical 
lightcurve usually consists two peaks (ZKM03 and ref- 
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erences there in), i.e., a forward shock peak which cor- 
responds to the epoch when the typical synchrotron fre- 
quency crosses the band (Sari et al. 1998; Kobayashi & 
Zhang 2003a) , and a reverse shock peak at which the flux 
achieves the maximum and starts to decay thereafter. 
This corresponds to the time when no more shock heat- 
ing is available and the shell starts to cool adiabatically. 
The condition for a reverse shock to exist is expressed in 
equations (|31[) and l|43l) . and in this paper we focus on 
the situation when such a condition is satisfied. We will 
denote the reverse (forward) shock peak times and peak 
fluxes as t p<r (t p j) and F UjPjT (F ViP j), respectively. 

For the a — case, the first shock crossing time and the 
shell deceleration time coincide, so that t PjT — t x = tdec- 
This is still the case when a is larger as long as the shell 
is in the thick shell regime (I). For thin shell cases (II 
and III), however, this is no longer the case, and t x and 
t^c separate from each other (eas.|5B].|5S|). For easy 
discussion, hereafter we define the time tdec as the reverse 
shock peak time, and its corresponding afterglow flux as 
reverse shock peak flux, i.e. 

tp,r — tdec 
F^,p,r=F u (tdec) (65) 

For t > t p>r , the shell cools adiabatically, and a decay- 
ing lightcurve results 7 . The t < t p>r case is a little more 
complicated. For the thick shell case (I), since t x co- 
incides with t Ptr , it is a rising lightcurve due to shock 
crossing. For the thin shell cases (II and III), as demon- 
strated above, the shell is heated at first shock crossing 
and remains heated until it is significantly decelerated. 
This results in a broad reverse shock peak which starts 
at t x and ends at t p>r . Much more detailed studies are 
needed to reveal the physics during this stage, but to 
a first-order estimate, in this paper we assume that the 
heating level between t x and t Pir is roughly the same, so 
that the lightcurve shows a plateau during this period. 

To calculate the synchrotron radiation flux, one needs 
to quantify the comoving random Lorentz factor of the 
leptons. We still take the convention to assume that the 
shock accelerated leptons have a single power-law distri- 
bution with the indices pf and p r for the forward and 
the reverse shocks, respectively, and that they occupy a 
fraction e e j and e e>r of the total thermal energy in the 
forward shock region (which is e 2 ) and in the reverse 
shock region (which is 63), respectively. For the reverse 
shock region, the lepton density may be enriched by the 
presence of pairs generated in the prompt emission phase 
(e.g. Li et al. 2003b). The pair- multiplicity parameter 

y= (N b + N ± )/N b > 1 (66) 

may be of order unity or mildly large in the low-er regime 
(depending on the compactness of the region when the 
prompt gamma-rays are emitted, e.g. Kobayashi, Ryde 
& MacFadyen 2002; Meszaros et al. 2002), and could 
be very large in the high-cr regime (e.g. Zhang & 
Meszaros 2002). In this paper, we are mainly focusing 
on the novel features introduced by the er parameter, 
and will take y ~ 1 in the following calculations. The 
y-dependences are included in the expressions and their 
implication will be discussed in ^4.41 

7 When a is very large, additional heating for the shell may 
still happen if the post-shock-crossing energy transfer process time 
scale is short enough. 



The minimum comoving electron energy in the re- 
gion "£" (2 or 3) is j e ,m,i = (^i/uMei/riimpC 2 ) 
g(Pi)(rn p /m e ), where g(p) = (p - 2)/(p - 1) (assuming 
p > 2). Taking the values at the first shock crossing time 
t x , one gets 

7e,m,r(*x) _ / £ e ,r 9r\ , 734 C*x) ~ 1 ^ f 70 -1 
/, \ — I I Ja , s 1 ~ ><-eJa 2 V ' 

7e,m,/C*x) \ye e jg f J J2{t x )-l 7x 

(67) 

where we have defined 

Tie = , (68) 

used 734 ~ 74/72(x) (which is valid for both thick and 
thin shells), and replaced 74 and 72 (x) by 70 (which 
means the initial Lorentz factor) and 7 X (the fireball 
Lorentz factor at the shock crossing time), respectively 8 . 
This allows the same notation system as in our previous 
work (ZKM03). 

We are more interested in the behavior at the reverse 
shock peak time (i.e. the deceleration time), t p . r . For the 
thick shell case, this is simply t x . For the thin shell case, 
after the first shock crossing, the shell is kept heated 
by multi-crossing of successive shocks. To first order, 
we can take the approximation that the heating level in 
the ejecta during the time period between t x and t p r 
is approximately constant, i.e., 63/77.3 oc t . During the 
same period the random Lorentz factor in the shocked 
ISM region also remains constant (since 72 ~ 74 before 
deceleration), we therefore also have 

7e,m,r(^p,r) 7e,m,r(^x) n J TO —1 rcn\ 

v — Tin ~ ~7 — jm ~ nja ^ y ' (69) 

le,m,f\ b p,r) !e.,m,f\'-X I l x 

which is valid for both the thick and thin shell cases. 

Let us denote the total electron numbers in the forward 
and reverse shocked region as N e j and N^ r , respectively. 
In the reverse shock region, the total lepton number (now 
includes pairs) is N e ^ r = N b + N± = yN b , where the 
definition in eq. lj6T)|> is used. According to ea. lBU|) . at 
the deceleration radius, the total energy in the forward 
shock region is defined by Ek alone. Although the bulk 
of the Poynting energy is expected to be transferred to 
the ISM eventually, shortly after the shock crossing and 
near the forward shock peak, this correction may not 
be significant. Below we will ignore this process in our 
calculations of the forward shock emission, but keeping 
in mind that the real forward shock emission level would 
increase with time, and may be much higher than our 
predicted level at later times. When we focus on early 
afterglow lightcurves, our calculations should be close to 
the real emission level (see more discussions in Zhang 
& Kobayashi 2005). A more careful treatment will be 
presented in a future work. 

In our approximated treatment, one can write 
Er = 7oc 2 (N b m p + N±m e ) — "f N b m p c 2 ~ 
Ne,fm p c 2 ['y(tdec)] 2 ~ Nejm p c 2 j x , where we have as- 
sumed y -C m p /m e , so that the total pair mass N±m e 

8 Strictly speaking, the last factor 70 /7J i n ea. 1671 should be 
(70 ~ 7x)/7x ■ The current approximation is valid as long as the 
reverse shock is mildly rclativistic, say, 734 > 1.5. For an even 
smaller 734 (which could be possible when a is large enough), there 
sho uld b e an add itional correction factor (less than unity) in both 
eqs. 1671 and 1691 . 
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is much smaller than the total baryon mass N m p . This 
gives 

-^e,r(^p,r) 7x 
N e j(tp, r ) 70 



(70) 



The characteristic synchrotron emission frequency is 
v m cx 7B7 2 m , the cooling frequency is 2/ c oc r y~ 1 B~ 3 t~ 2 , 
and the peak specific flux is cx jBN e , where 7 is the 
bulk Lorentz factor. For the thin shell case, we also make 
another approximation that B3 keeps constant from t x 
to ^ec (so that TZs{tx) — ^-s^dec) — T^-b)- Similar to 
Kobayashi & Zhang (2003a) and ZKM03, we can finally 
derive the following relations at t ( i ec - 

^m,r(^p,r) 



^c,r(^p,r) 
^c,f (^p,r) 
-fl^,m,r(^p,r) 



where 



_7x . f 7c 2 
7 = — = min 7o, — 
7o V 7o 



< 



7c 



(71) 
(72) 
(73) 

(74) 



Although there are in principle many cases of the 
reverse shock emission lightcurves (Kobayashi 2000), 
within the reasonable parameter regime the lightcurve 
behavior only has two variations depending on whether 
1Z V > 1 or H v < 1 (ZKM03) where 

n v = J . (75) 

In both cases, the ratio between the two peak-time fluxes 

F v . 



K F = 



(76) 



(77) 



and the ratio between the two peak times 

t 

can be expressed in terms of 7, TZb and IZv, respec- 
tively, for the cr = case (ZKM03). Here the forward 
shock peak time t p j corresponds to the epoch when v m j 
crosses the band 9 . 

Below we repeat this process, but focus more on cor- 
rections introduced by the a factor. To further simplify 
the problem, we first estimate the numerical value of TZ^. 
Due to the complication introduced by the a parameter, 
one can not coast 1Z V into a simple expression as in the 
a = case (e.g. eq.[24] in ZKM03). In any case, us- 
ing eq.J7TJ and the standard expression for v m ,f(t) (e.g. 
eq.[l] in Kobayashi & Zhang 2003a), one can write 



_ 2 ( kE K 



\ 10 52 erg, 

/ € bj y 1 ' 2 f^iy 2 ( 9j_ s ~~ 
Vo.001/ \0.lJ U/3, 

/ tde 



-1/2 / . x 2 

' 72 \ 

100/ 



3/2 



V100 s 



l + z 



-1/2 



(78) 



9 When a is large enough, additional correction factor to t p j is 
needed, but this factor is small if the energy injection time scale is 
long enough. 



For 7 < 7c ~ 125 (eq.[H]), TZ e ~ 1, y > 1, and /„ < 
1, the above equation therefore essentially always gives 
1Z V > 1. In the following discussions, we will not discuss 
the 1Z V < 1 case any further (which was also discussed 
in ZKM03). 

The reverse shock emission lightcurve in the 1Z V > 1 
case is simple. The lightcurve initially rises and reaches 
the peak at t x . The flux level then keeps essentially 
constant until t v . r (for the thick shell case, both time 
scales coincident, so that there is no broadened peak), 
and starts to decay after t p ^ r . The temporal indices of 
each segment of the lightcurve is also well-defined. In 
the rising part of the lightcurves, since all the correc- 
tion factors introduced by the a parameter are essentially 
time- independent, the corrections essentially do not in- 
troduce extra time-dependence on the typical frequencies 
and the peak flux of the synchrotron radiation in the re- 
verse shock. The rising lightcurves essentially remain 
unchanged as the a = has been derived by 

Kobayashi (2000). This gives a ^ 1/2 temporal index 
for the thick shell case, and a ~ 5 temporal index for 
the thin shell case 10 . After the deceleration time, the 
shell cools. The optical band is typically in the regime of 
v m ,r{tx) < V R < v c , r (t x ). After the deceleration time, 
one has v m ^ r cx £~ 3 / 2 , F u>mtT cx i -1 (Kobayashi 2000) 11 . 
Thus the temporal decay index (i.e. F v cx t~ a ) is 



3p r 



1 



(79) 



where p r is the electron power-law index in the reverse 
shock region. 

For t > t p<r , one has v m j cx i -3 / 2 , F v<m j cx t° 



(Meszaros & Rees 1997a) and v„ 



+-3/2 F 



t^ 1 (Kobayashi 2000). Using the definitions of TZ V , TZf 
and lZ t (eas. |78lf7| ) as well as eqs. (|7T|l and (J75J, one 

can derive 13 

K t = 7 4/3 7^ 2/ X~ 2/ W 7 V /3 / a - 4/3 ) , (80) 

TZ F ^jn B TZ- 2 ^/ 3 (y) . (81) 

These are valid for all the three parameter regions in 
Fig|U Comparing with eqs. (12, 13) in ZKM03, the extra 
correction factors are presented in parenthesis. Notice 
that the correction factors TZ e and y should also exist 
in a = case, but we have previously assumed them to 
be unity. The extra c-dependent correction factors are 

f a 4 ^ 3 and 1Z U (which is modified by the a parameter 
through many factors, e.g. TZb, fa, Ek and tdec, see 
eq.[?E!). 

4.4. Sample lightcurves 

We now calculate the typical early optical afterglow 
lightcurves for various parameter regimes. Equation l|40|l 

10 Detailed numerical calculations result in non-power law be- 
havior in the rising lightcurve (Fan et al. 2004a). 

11 A more detailed discussion such as that presented in 
Kobayashi & Sari (2000) leads to the similar conclusion. 

12 Notice again that here we have assumed that the energy trans- 
fer time scale from a Poynting flux to the kinetic energy of the ISM 
long enough. This forward shock emission level should be regarded 
as a lower limit when the energy transfer process is taken into 
account. 

13 In ZKM03, we have defined TZp and IZt at i x , but in a = 
case one has tx = idee- For the case of an arbitrary a, the 
deceleration time t^ ec is more fundamental to define the problem. 
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states that the initial afterglow energy, which is essen- 
tially the kinetic part of the total energy, decreases with 
a given a constant total energy E = Ex+Ep. At high-<7, 
not only the reverse shock flux level drops, the forward 
shock flux level shortly after the shock crossing also de- 
creases steadily. At later times, the forward shock level 
would increase due to magnetic energy injection. Since 
we are focusing on the early afterglow emission, this ef- 
fect will be neglected in the following discussions. To 
explore the c-effect, we fix the total energy of the fire- 
ball so that Ek decreases with increasing a. To simplify 
the calculations, we assume lZ e ~ 1 and y ~ 1. The 
input parameters we adopt include E^ 2 = 1, 7o = 150, 
n = 1, e eJ = 0.1, e B j = 0.001, p f = 2.2 and z = 1 
(with the standard cosmological parameters Q\ ~ 0.7, 
fl m ~ 0.3 and H ~ 70 km s^ 1 Mpc -1 ). This gives the 
forward shock peak time and flux (Sari, Piran & Narayan 
1998; Kobayashi k, Zhang 2003a) 

t pJ ~ 1000 s (82) 
1.7(1+ tr)- 1 mJy 

[m fl ~ 15.6 + 2.5 log(l + cr)] . (83) 

We also have t 7 = [(3E /A-K^nmpC 2 ) 1 / 3 /2^c](l + z) ~ 
60 s. (For this calculation, we have added in all the 
precise coefficients previously neglected.) We then take 
two typical values of GRB durations. For the first case, 
we take T = 100 s. When a is small the burst is in the 
thick shell regime (region I). As a increases, the burst 
is in the non-spreading thin shell regime (region II). For 
the second case, we take T = 20 s. The burst is always in 
the thin shell regime for any a value, but transform from 
the spreading thin shell regime (region III) to the non- 
spreading thin shell regime (region II) when a is large 
enough. 

For each T value, we calculate both the reverse shock 
and the forward shock lightcurve for several values of 
cr, i.e. a = 0,0.001,0.01,0.1,1,10,100 (FigEJ), as long 
as the condition for the existence of the reverse shock 
fea.|43|l is satisfied. The procedure of our calculation is 
the following. First, with t 7 , T and the assumed cr, one 
can judge which parameter region the burst is in. With 
this information one can then calculate t^ ec = t p>r and 
7Zt for both the thick and thin shell regimes, as well as 
t x for the thin shell case. Next, we calculate TZb with 
the assumed a value (FigEJ. F° r the thick shell case, 
we use the Zb value for 734 ~ 1000 since Zb is insensi- 
tive to 734 when it is large. For the thin shell case we 
use the Zb value for 734 ~ 1.5, exclusively . One can 
then solve 1Z U (eqJHOj), and then use the value of 1Z V to 
calculate TZf (eq.jSU), and hence F ViP)T . Since we know 
the temporal indices of the reverse shock lightcurve dur- 
ing the rising (~ 1/2 for thick shell and ~ 5 for thin 
shell, Kobayashi 2000) and the decaying phase (~ —1.9 
for p r = 2.2), the reverse shock lightcurve can be calcu- 
lated once tp^r = tdec and F UtPtT are known. For the thin 
shell case, with the current approximation, we roughly 
keep F v a constant between t x to td ec , both of which are 
known. For the forward shock emission, the temporal 
index is 3(1 — p/)/4 (~ 0.9 for pj = 2.2) after the peak 

14 Notice that in reality, when a is very large, 734 could be much 
closer to unity. In such cases, the reverse shock peak flux should 
be further suppressed. 
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Fig. 7. — Sample early afterglow lightcurves for GRBs with an 
arbitrary magnetization parameter cr. Following parameters are 
adopted. £52 = 1, 70 = 150, n = 1, e e j = 0.1, e B j = 0.001, pj = 
2.2, and z = 1. We also take lZ e ~ 1 and y ~ 1. Both the forward 
shock and the reverse shock emission components are calculated 
and they are superposed to get the final lightcurve. For the forward 
shock emission component we assumed that the time scale for the 
energy transfer from the Poynting energy to the afterglow energy is 
long enough, so that the forward shock level is defined by Ex only, 
and its level decreases with a. This approximation is good shortly 
after the reverse shock crossing. At later times, the real level could 
be progressively higher than this level and the calculation should be 
regarded as a lower limit. Lightcurves are calculated for different 
cr values. Thick solid: a = 0; thin dash-dotted: cr = 0.001; thin 
dashed: cr = 0.01; thin dotted: cr = 0.1; thin solid: cr = 1; thick 
dash-dotted: cr = 10; thick dashed: cr = 100. 734 = 1.5 has been 
assumed for thin shell regimes. For high-cr cases, 734 is closer to 
unity, and the reverse shock peak flux sho uld be further suppressed, 
(a) T = 100 s case. According to eg. 1431 . the reverse shock exists 
when cr < 200. (b) T = 20 s case. The condition for the existence 
of the reverse shock is cr < 20. 



time, and is 1/2 before the peak time (but after the de- 
celeration time) 15 . Given F p j (which is dependent on a 
fea.|83|l. the forward shock lightcurve is also calculated. 

Some sample R-band early afterglow lightcurves are 
presented in Fig0 with the contributions from both the 
reverse and the forward shocks superposed. For the for- 
ward shock emission component we assume that the time 
scale for the energy transfer from the Poynting energy to 
the afterglow energy is long enough, so that the forward 
shock level is defined by Ek only, and its level decreases 
with cr. This approximation is good shortly after the re- 
verse shock crossing. At later times, the real level could 
be progressively higher than this level and the calcula- 

15 Before the deceleration time, the forward shock lightcurve 
should have different temporal slopes. During the shock crossing, 
we have 72 oc t° for thin shells and 72 on t -1 / 4 for thick shells. 
Using the standard synchrotron radiation analysis (e.g. Sari et 
al. 1998), the forward shock emission temporal slope is 3 and 4/3 
for the thin and thick shell cases, respectively. Between the shock 
crossing time and the deceleration time in the high-cr thin shell 
case, the temporal slope is fiat. 
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tion should be regarded as a lower limit (see Zhang & 
Kobayashi 2005 for more explanations). In FigEJi, the 
cases for T — 100 s are calculated. According to ea. l|43|) . 
a reverse shock exists as long as a < 200. We therefore 
calculate the lightcurves up to a = 100. We can see that 
for (j < 1, the parameters are in the thick shell (relativis- 
tic reverse shock) regime. When a increases from below, 
the contrast between the reverse and forward shock peak 
fluxes (i.e. TZp) increases steadily. Since the forward 
shock emission level does not change much when a < 1, 
the reverse shock peak flux increases steadily with a. At 
even higher a values, the reverse shock peak flux drops 
steadily with cr. In the mean time, the burst enters the 
thin shell regime so that the separation between t x and 
tdec becomes wider, and the reverse shock emission has a 
broader peak. In Fig 03, the cases for T = 20 s are calcu- 
lated. According to ea.l|43|). a reverse shock exists as long 
as a < 20, and we calculate the lightcurves up to a = 10. 
The shell is in the thin shell regime throughout the whole 
a range calculated. The transition from spreading to 
non-spreading thin shell regime does not bring any no- 
ticeable signature in the lightcurves. Again, the reverse 
shock peak flux increases with a initially (when a < 0.1), 
and starts to decrease when a > 0.1. The reverse shock 
peak is broad, but the separation gradually shrinks due 
to the decrease of the (Q -1 ^ 2 — Q~ x ) parameter (eq.jSSJ). 
Throughout our calculations, TZ V remains larger than 25 
(up to ~ 1000 for a = in the T = 100 s case), so that 
our treatment by neglecting TZ V < 1 regime is justified. 

We notice several interesting features from our results. 
First, the reverse shock component is still noticeable even 
with <J>1 (until reaching several tens or even hundreds 
when the condition 1|43[1 is no longer satisfied) . The abso- 
lute reverse shock peak flux increases with a initially, but 
drops steadily when a > 1. Second, the forward shock 
emission level right after shock crossing also drops with 
cr. This is because only the kinetic energy of the bary- 
onic component (Ek) defines the afterglow level after the 
shock crossing time. The forward shock level will increase 
later due to the transfer of the remaining magnetic en- 
ergy into the medium. One then expects an initially dim 
early afterglow for a high-cr flow, which would be bright- 
ened at later times. If GRB prompt emission is due to 
magnetic dissipation (e.g. Drenkhahn & Spruit 2002), 
and if a is still high in the afterglow phase (e.g. ~ 10), 
one may account for the very large apparent GRB effi- 
ciencies inferred from some GRBs (e.g. Lloyd-Ronning 
& Zhang 2004). Such a picture may be also relevant 
to the recent December 27 giant flare afterglow from 
the soft gamma-ray repeater 1806-20, for which a very 
high gamma-ray efficiency is inferred (Wang et al. 2005 
and references therein). Third, the broad reverse shock 
peak is a novel feature identified in the high-er model, it 
can be used to diagnose the existence of a Poynting-flux- 
dominated flow. The physical origin of the broad peak 
is that a high-cr value leads the decoupling of the shock 
crossing radius i? x and the deceleration radius i? 7 , so 
that multi-crossing of a series of successive shocks leads 
to continuous heating of the cjecta shell before cooling 
starts. 

In the above calculations, y = 1 has been adopted 
(i.e. we assume that the pair fraction is negligible in the 
ejecta). In some cases, especially in the high-cr regime, 



y could be much larger than unity. It would be essen- 
tial to investigate the y-dependence of the current analy- 
sis. Solve TZ V from eci. (|SU)l and submit it to TZf, we find 
TZf oc y( 7 ~ 4a )/ 3 , which is oc y' 1 ^ for a = 2. We can see 
that a larger y will lower the reverse-to-forward shock 
peak flux contrast, although the dependence is mild. The 

TZf factor is more sensitive to lZ e (i.e. oc TZ^ a ), 
but assuming a similar shock acceleration mechanism, 
TZ e may not deviate too much from unity. 

Our results can be directly compared with the early 
afterglow data of the four bursts whose such information 
is available so far. The case of GRB 990123 (Akerlof et 
al. 1999) is consistent with a flow with 0.1 < a < 1 in 
which regime TZf is large and the reverse shock peak is 
not broadened. The observed bright afterglow also ar- 
gues against a higher-cr flow. The case of GRB 021211 
(Fox et al. 2002b; Li et al. 2002a) also shows a large 
TZf i which also suggests that a > 0.1. For GRB 021004, 
Kobayashi & Zhang (2003a) attempted to fit the data 
with the a = 0, TZb = 1 model. Another data point 
at an earlier epoch after the burst trigger reported by 
Fox et al. (2003a) makes that model more difficult to 
fit, and it has been attributed to a continuously energy 
injection (Fox et al. 2003a) or to the emission from a 
wind- type medium (Li & Chevalier 2003). However, us- 
ing the theory developed in this paper, the data may be 
consistent with a high-cr flow (e.g. cr^lO), so that the ex- 
tended early afterglow emission could be interpreted as 
the combination of the broad reverse shock peak and the 
gradual transfer of the Poynting energy into the after- 
glow energy. In the high-cr regime, the (734 — 1) is quite 
small during the shock crossing, which will lower the re- 
verse shock peak flux and TZf- Similarly, a broad early 
afterglow bump was identified in GRB 030418 (Rykoff et 
al. 2004) which challenges the conventional reverse shock 
model but may be consistent with a high-cr flow (weak 
or no reverse shock component) . Although detailed mod- 
eling is needed (we plan to do it in a future work), we 
tentatively conclude that all the current early optical af- 
terglow data may be understood within the theoretical 
framework developed in this paper, if a is allowed to vary 
for different GRB fireballs. 

5. CONCLUSIONS AND DISCUSSION 

We have derived a rigorous analytical solution for the 
rela tivist ic 90° shocks under the ideal MHD condition 
fca. |A16j L Generally, the solution depends both on the 
magnetization a parameter and the Lorentz factor of 
the shock, 712. The solution can be reduced to the 
Blandford-McKee hydrodynamical solution when a = 0, 
and to the Kennel-Coroniti solution (which depends on 
cr only) when the 721 S> 1. Our generalized solution 
can be used to treat the more general cases, e.g. when 
the reverse shock upstream and downstream are mildly 
relativistic with each other. Since GRBs invoke a shell 
with finite width, this latter possibility is common (e.g. 
the parameter space for thin shell greatly increases in 
the high-cr regime), so that our generalized solution is 
essential to deal with the GRB reverse shock problem. 

Several interesting conclusions emerge from our analy- 
sis. (1) Strong shocks still exist in the high-cr regime, as 
long as the shock is relativistic. Figs|2:,d indicate that as 
cr increases, both the downstream "temperature" eijni 
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and the "shock compression factor" nijri\ decreases with 
respect to the cr = values. However, the suppression 
factors in both cases are only mild (a factor of ~ 0.5), 
and they saturate when er 3> 1. In the relativistic shock 
regime, the results are actually consistent with Kennel & 
Coroniti (1984). However, these authors did not calcu- 
late the suppression factor with respect to the a = case, 
and did not explore further into the high-cr regime, so 
that their results leave the impression to the readers that 
the shock is completely suppressed when a reaches higher 
values. For typical GRB parameters, we found that a re- 
verse shock still exists when cr is as high as several tens or 
even hundreds. When the reverse shock exists, its emis- 
sion level decreases when cr gets higher. This is not only 
because the reverse shock becomes weaker since 734 gets 
close to unity in the high-cr regime, but also because the 
total kinetic energy in the flow (which is the energy reser- 
voir for shock dissipation) gets smaller given a same total 
energy. (2) When discussing ejecta-medium interaction, 
somewhat surprisingly, some important parameters, such 
as F and Ca, are very insensitive to the reverse shock 
Lorentz factor, 734, and can be regarded as the function 
of a only f H2.3l M3.3[l . This greatly simplify the problem, 
and is essential to characterize the parameter regimes. 
(3) The triple coincidence of the first three critical lines 
in ea. Q49[) is very crucial for a self-consistent description 
of the problem (JOJ. 

Comparing with the conventional hydrodynamical 
treatment, we reveal several novel features for the early 
lightcurves. First, as a increases, the reverse shock peak 
flux increases rapidly initially, reaching a peak around 
a ~ (0.1 — 1), and starts to decrease when a > 1. Second, 
due to the inability of tapping the Poynting flux energy 
during the shock crossing process, the fireball decelera- 
tion radius for the thin shell case decreases as a increases 
[oc (1 + cr)~ 1//3 ]. The forward shock emission level is also 
lower right after shock crossing. Third, in the high-cr thin 
shell regime, the reverse shock peak is broadened due to 
the separation of the shock crossing radius and the de- 
celeration radius. This is a signature for a high-cr flow, 
which can be used to diagnose the magnetic content of 
the fireball. Fourth, as a becomes large enough (larger 
than several tens or several hundreds) , the condition for 
forming a reverse shock is no longer satisfied, and there 
should be no reverse shock component in the early af- 
terglow lightcurves. This could be consistent with very 
early tight optical upper limits for some GRBs, such as 
the recent Swift dark burst GRB 050319a (Roming et al. 
2005). In summary, the above new features allows the 
current theory to potentially interpret known GRB early 
afterglow cases collected so far as well as the case of the 
dark bursts, if one allows a to vary in a wide enough 
range (say, from 0.01 to 100). The Swift GRB mission, 
launched on November 20, 2004, is expected to detect 



many early optical afterglow lightcurves with the UV- 
optical telescope on board. We expect to further test our 
theoretical predictions against the abundant Swift data, 
and to systematically diagnose the magnetic content of 
GRB fireballs. 

If the GRB ejecta is indeed magnetized, as inferred 
from the early afterglow data, the internal shocks should 
also be corrected by the magnetic suppression factor. 
This aspect has been investigated by Fan, Wei & Zhang 
(2004b) recently. 

Throughout the paper, we have treated the problem 
under the ideal MHD limit. In the high-cr case, strong 
magnetic dissipation may occur. The magnetic dissipa- 
tion effect has been included in the internal shock study 
of Fan et al. (2004b). Our treatment in this paper 
presents a first order picture to the early afterglow prob- 
lem (see also Fan et al. 2004a, whose treatment in the 
mildly magnetized regime is consistent with ours), and 
further considerations are needed to fully delineate the 
physics involved. Also our whole discussion is relevant 
when a reverse shock is present. It does not apply to 
the regime for an even higher a value (e.g. Lyutikov & 
Blandford 2003). Finally, as discussed in §4.1, further 
studies on particle acceleration in MHD shocks are es- 
sential to give a more accurate calculation on the reverse 
shock emission in the high-cr regime. 

We only discussed one type of the medium, i.e., assum- 
ing a constant medium density, typically for the interstel- 
lar medium. In principle, the medium density can vary 
with distance from the central engine. In particular, a 
wind-type medium, characterized by the n cx R~ 2 profile, 
has been widely discussed. Our MHD shock theory could 
be straightforwardly used to the wind case to study the 
reverse shock emission in combination with the previous 
pure hydrodynamical treatments (Chevalier & Li 2000; 
Wu et al. 2003; Kobayashi & Zhang 2003b; Kobayashi, 
Meszaros & Zhang 2004, see Fan et al. 2004a for a pre- 
liminary treatment). 
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APPENDIX 

DERIVATION OF THE SOLUTION OF THE RELATIVISTIC 90° SHOCK JUMP CONDITIONS 

Lorentz transformations 

Given the definitions of 7^ , /3ij and Uij, the following Lorentz transformations are frequently used in the derivations. 



2s 1 - W21 



(Al) 
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a 02s + 021 

Pis — 



02V 

12s - 
lls- 



1 + 02s021 
Pis ~ 02s 

' 1 - 01S02S 
: 7ls72l(l - 01s02l) 
: 72s72l(l + 02s02l) 
721 =7ls72s(l - 0U02s) 
W2s=7ls72l(/3ls ~ 02l) 
Uls=j2sj2l(02s + 02l) 
"21 =7ls72s(/3ls - 02s) 

R R - U21 

01s — 02s — 

7ls72s 

Derivation of equation 



Let us define (Kennel & Coroniti 1984) 



Y 



B2s 

Bis 



l2sUls 
llsU 2s 



0U 
02s'' 



the equations © and Q) can be re-written as 

7i s /Ul[l + (1 - Y)a] =j2sti2 

(7 



u ls ni 



1 



2/3?.. 



= U 2 s^2 



P2 



n 2 U2s 



Multiply ca. (|A12l) by 0i s and substitute the resultant formula into ea. l|A13|) . one can derive 



n 2 m p c 



e 2 



»2 



(71,7241 + (1 - Y)a\ - Ui s U 2 s[l + \{0i^ - 02s 2 )]) Ml 



721 



'-21 



a iii. 



2ui s u 2s 

The equation (JSJ) in the text can be then derived straightforwardly. 



Solving u 2 



2s 



Combining ea. (|A12|l and the definition of [12 (eq.[5]), one can derive 



^[1 + (1-F)a] = l + f(72i 

72s 



1) 



2u ls U 2s 



This turns out to be a three-order equation of x 



l 2s> 



i.e. 



where 



Ax 3 + Bx 2 + Cx + D = 0, 

A = f(2-f)( 721 -l) + 2 

B = -(721 + 1) [(2 - f)(f 72 2 1 + 1) + f (f - 1) 72 

-(721 - 1) [f (2 - f )( 7 f 1 - 2) + (272! + 3) 



C = ( 7 2i + 1) 



r(i- I )(7 2 2 i-i) + i 



+(721 
+(721 



1) [2721 - (2-r)(r 72 i - 1) 

l)(72i-l) 2 (f-l) 2 



^ = -(72i-l)(72i + l) 2 (2-f) 2 T . 

For r = 4/3, the four coefficients could be written equivalently as 
A = 8721 + 10 

B = -(721 + l)(8 7 |i + 4 7 2i + 6)cr - (721 - 1)(8 7 | 1 + I8721 + 11) 

C= (721 + l)(8 72 2 i + l)^ 2 + (7 2 2 i ~ l)(1072i + 6)ct + (721 + l)(72i - I) 2 

^ = -(72i-l)(72i + l) 2 ^ 2 • 
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When a = 0, the equation l|A16(l is reduced to 



17 



[* - (ill ~ 1)] {[f (2 - f )( 721 - 1) + 2]x - (721 - l)(f - l) 2 } 



(A25) 



which gives equation l|12l) besides the other two non-physical solutions U2 S — and w 2s = u 2 \. 



721 ^ 1 limit 

When 721 ^> 1, the x 3 term is a small quantity and is negligible. The equation (|A16|) is reduced to 



f(2-f)0 + l)a; 2 - 



f (1 - -)(7 2 + (f 2 - 2f + 2)cr + (f - l) 2 



x + (2-r) 2 - = o 



(A26) 



This gives the solution H17|) in the text (when the non-physical solution is neglected). 

NOTATION LIST 

subscript 1 upstream (§2 and Appendix A), or unshocked medium (§3 and §4) 
subscript 2 downstream (§2 and Appendix A), or shocked medium (§3 and §4) 
subscript 3 shocked ejecta 
subscript 4 unshocked ejecta 
shock 

speed of light 

internal energy density in region i(— 1, 2, 3, 4) 



subscript 
c 
e, 

fa 
fb 
fc 

9f 
9r 
I 

m e 

m p 

m 

Pi 

Pb,i 

Pf 

Pr 
tdec 

l vJ 

tp r 

tx 
Ujj 

X 

y 

z 



correction factor of (eijniTdyC ) normalized to the a — value 
correction factor of (712/ni) normalized to the a = value 
magnetic-to-thermal pressure ratio plus 1 
(p f - 2)/(p f - 1) 

( Pr - 2)/{p r - 1) 

Sedov length 
electron rest mass 
proton rest mass 

baryon number density in region i{= 1, 2, 3,4) 
thermal pressure in region i(= 1,2,3, 4) 
magnetic pressure in region i(= 1, 2, 3, 4) 
electron power law index in the forward shock 
electron power law index in the reverse shock 
deceleration time measured by the observer 
J R 7 /C 7 7|c (eq.pj) 

emission peak time of the forward shock component 

emission peak time of the reverse shock component 

shock crossing time measured by the observer 

four speed in the region i{= 1, 2, 3, 4) in the rest frame of j(— 1, 2, 

pair multiplicity parameter 
redshift 



3,4,s) 
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A, B, C, 


D coefficients to solve the equation for u\ s 




B, 


comoving magnetic field in the region i(— 1, 2, 3, 4) 




Bi S 


magnetic field in the region i(— 1,2,3, 4) in the rest frame of the shock 


T~> 

B f 


■ j*P11*j1P 1111 • 

comoving magnetic field m tne forward shocked region 




77) 

D T 


■ j * f* 1 1 * j 1 ill 

comoving magnetic field m tne reverse shocked region 




Oat 


j ■ f 1 1 Tt • 1 1 1 1 1 1 r\ 

correction lactor to with respect to the a — case 






correction factor to Ka with respect to the a = case 




/ 1 


correction factor to i? 7 with respect to the a — case 




c 

t 


1 1 C 11*011 

shock frame electric field 




TP 


•i *jj1 PjI f~*l"ll 

isotropic total energy of the fireball 




771 


isotropic kinetic energy of the fireball 




TP 

Hip 


isotropic Foyntmg-liux energy of the hreball 




TP 

t 


the product of / a , f b and f c 




771 


specific flux at the frequency v 




TP 

" v,m,f 


• ij ■ ■ •pn'jip ill 

maximum synchrotron emission specific flux m the forward shock 






maximum synchrotron emission specific flux in the reverse shock 






1 r% p . ■! p 111 j • j • / 

peak flux for the forward shock emission component m certain (e.g 


;. R) band 


TP 


peak flux for the reverse shock emission component in certain (e.g. 


R) band 


TT 


Hubble constant 




a 7 

M 


mass in the ejecta 




-("ISM 


mass of the interstellar medium collected by the shock 




A 7" 

iV& 


total baryon number in the shell 




A T 

N± 


total electron-positron pair number in the shell 




AT 

N eJ 


lepton (electron) number in the forward shock 




AT 

N e ,r 


lepton (electron and pairs) number in the reverse shock 




Q 


a parameter introduced to categorize the parameter regimes (defined by eq.joOj 


T~) 

E 


radius from the central engine 




TIb 


reverse-to-forward comoving magnetic field ratio 






reverse-to-forward ratio of the e e g parameter 




Ti F 


reverse-to-forward peak flux ratio 




Tli 


forward-to-reverse peak time ratio 






the ratio between vr and v m ,r{tdec) 




77) 

Kn 


radius where the reverse shock becomes relativistic 






radius where the ejecta shell starts to spread 






radius where the reverse shock crosses the ejecta shell 




77) 

JX>y 


radius where the fireball collects I/70 rest mass of the fireball 




T 


central engine activity time scale 




U B ,o 


initial comoving magnetic energy 




U b 


comoving magnetic energy after shock crossing 




X 


an intermediate parameter introduced in ea.(|17|) 




Y 


ratio between B^ s and B\ s 
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r adiabatic index 

A shell width in the lab frame 

Ao initial shell width in the lab frame 

Q m cosmology mass density parameter 

J7a cosmology A density parameter 

a temporal decay index of the reverse shock emission component after peak time 

0ij dimensionless velocity of region i(= 1, 2, 3, 4) in the rest frame of j(= 1, 2, 3, 4, s) 

7 an equivalent Lorentz factor defined in ea. l74|) 

7y Lorentz factor of region i(— 1, 2, 3, 4) in the rest frame of j(= 1, 2, 3, 4, s) 

7i Lorentz factor of region i(= 2,3,4) in the rest frame of the circumburst medium 

70 initial Lorentz factor of the fireball, 70 = 74 

7c critical initial Lorentz factor that separates thick vs. thin shell regimes feci. |52j) 

7x hreball Lorentz factor at the shock crossing time 

7e,m,/ electron minimum Lorentz factor in the forward shock 

7e,m,r electron minimum Lorentz factor in the reverse shock 

e e j electron energy equipartition parameter in the forward shock 

e e ,r electron energy equipartition parameter in the reverse shock 

€bj magnetic energy equipartition parameter in the forward shock 

£B,r equivalent magnetic energy equipartition parameter in the reverse shock 

fii specific enthalpy in region i(= 1, 2, 3, 4) 

v c .f forward shock synchrotron cooling frequency 

i/ Cjr reverse shock synchrotron cooling frequency 

v m.j forward shock synchrotron typical frequency 

v m ,r reverse shock synchrotron typical frequency 

vji R-band frequency 

£ a parameter defined in ea. H45(l 

£0 the £ value when A = Ao 

a magnetization parameter as defined in eq.JTJ) 
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